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Abstract 

We propose a method for determining 7 using — > DK^ decays followed by a multibody D 
decay, such as D — > K$ 7r~7r + , D — ► Kg K~K + and D — ► 7r _ 7r + 7r°. The main advantages of the 
method is that it uses only Cabibbo allowed D decays, and that large strong phases are expected 
due to the presence of resonances. Since no knowledge about the resonance structure is needed, 7 
can be extracted without any hadronic uncertainty. 
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I. INTRODUCTION 



The theoretically cleanest way of determining the angle 

1 = ^g{-v ud v: b /v cd v: h ) ) (i) 

is to utilize the interference between the b — > cus and b — > ucs decay amplitudes [1-12]. Be- 
cause these transitions involve only distinct quark flavors, there are no penguin contributions 
to these decays. In the original idea by Gronau and Wyler (GW) [3] the B ± — > D C pK ± 
decay modes are used, where D C p represents a D meson which decays into a CP eigen- 
state. The dependence on 7 arises from the interference between the B ± — > D°K ± and 
— > D If ± decay amplitudes. The main advantage of the GW method is that, in principle, 
the hadronic parameters can be cleanly extracted from data, by measuring the B ± — > D°K ± 
and B ± — > D K ± decay rates. 

In practice, however, measuring 7 in this way is not an easy task. Due to the values of the 
CKM coefficients and color suppression, the ratio between the two interfering amplitudes, 
r B [see Eq. (4)], is expected to be small, of order 10% — 20%. This reduces the sensitivity to 
7, which is roughly proportional to the magnitude of the smaller amplitude. In addition, if 
the strong phases vanish, measuring 7 makes use of terms of order r\. In contrast, if a large 
strong phase is involved in the interference, there is a sensitivity to 7 at order r B with most 
methods. Thus, in general, having large interfering amplitudes with large relative strong 
phases is a favorable situation. 

Since the hadronic parameters are not yet known, it is still not clear which of the proposed 
methods is more sensitive. In addition, all the methods are expected to be statistically 
limited. It is therefore important to make use of all modes and methods, as well as to 
try to find new methods. Any new method that is based on "unused" decay channels 
increases the total statistics. Moreover, many of the analyses are sensitive to common 
hadronic parameters, for example, r B . Combining them will increase the sensitivity of the 
measurement by more than just the increase in statistics. 

Here we study the possibility to use B ± — > DK ± , followed by a multibody D decay, in 
order to cleanly determine 7. While this idea was already discussed in [5], most of our results 
and applications are new. For the sake of concreteness, we concentrate on the D — > K s 
decay mode. The advantage of using such decay chains is threefold. First, one expects large 
strong phases due to the presence of resonances. Second, only Cabibbo allowed D decay 
modes are needed. Third, the final state involves only charged particles, which have a higher 
reconstruction efficiency and lower background than neutrals. The price one has to pay is 
that a Dalitz plot analysis of the data is needed. We describe how to do the Dalitz plot 
analysis in a model-independent way, and explore the advantages gained by introducing 
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verifiable model-dependence. The final balance between the advantages and disadvantages 
depends on yet-to-be-determined hadronic parameters and experimental considerations. 



II. MODEL INDEPENDENT DETERMINATION OF 7 

As we shall show in this section, to perform a model independent determination of 
the angle 7 one needs to measure the two CP-conjugate decay modes, B ± — > DK ± — > 
(K S 7r~TT + ) D K ± and to perform a Dalitz plot analysis of the K s tt~it + final state originating 
from the intermediate D meson. (In the following discussion we neglect D° — D° mixing, 
which is a good approximation in the context of the Standard Model. See appendix A for 
details.) 

Let us first focus on the following cascade decay 



B~ -> DK~ -> (K s tt-tt + ) d K-, (2) 



and define the amplitudes 



A(B~ ^D°K-) = A B , (3) 
A(B~ -> D°K~) = A B r B e t{6B -^. (4) 

The same definitions apply to the amplitudes for the CP conjugate cascade B + — > DK + — > 
{Kg n + n^) B K + , with the change of weak phase sign 7 — > —7 in (4). Since we have set the 
strong phase of A B to zero by convention, 5 B is the difference of strong phases between the 
two amplitudes. For the CKM elements, the usual convention of the weak phases has been 
used. (The deviation of the weak phase from —7 has been neglected, as it is suppressed 
by the factor A 4 ~ 2 x 1CT 3 , with A being the sine of the Cabibbo angle.) The value of 
\A B \ is known from the measurement of the B~ — > D°K~ decay width using flavor specific 
decays of D° and the precision of its determination is expected to further improve [13]. The 
amplitude A(B~ — > D°K~) is color suppressed and cannot be determined from experiment 
in this way [4]. The color suppression together with the experimental values of the ratio of 
the relevant CKM elements leads to the theoretical expectation r B ~ 0.1 — 0.2 (see recent 
discussion in [11]). 

For the three-body D meson decay we define 

A D (s 12 , s 13 ) = A 12>13 ee A(D° - K s ( Pl )n-(p 2 )n + (p 3 )) 

n (5) 

= A(D° ^K s ( Pl )7T + (p 2 )7r-(p 3 )), 

where s^- = (pi +pj) 2 , and Pi,p 2 ,p 3 are the momenta of the A" 5 ,7r~,7r + respectively. We 
also set the magnitude A 12jl3 > 0, such that #12,13 can vary between and 2n. In the last 



3 



equality the CP symmetry of the strong interaction together with the fact that the final 
state is a spin zero state has been used. With the above definitions, the amplitude for the 
cascade decay is 

A{B- - (K s n-n + ) D K-) = A B V D (A D (s 12 , s 13 ) + r B e^ SB ^A D (s 13 , s 12 )), (6) 

where Vd is the D meson propagator. Next, we write down the expression for the reduced 
partial decay width 

dt(B- - (K s n-n+) D K~) = (a? 2>13 + r\ A^ 12 

\ (7) 

+ 2r B TZe [A D {s 12 , s 13 ) A* D (s 13 , s 12 ) e' 1 ^^] jdp, 

where dp denotes the phase space variables, and we used the extremely accurate narrow 
width approximation for the D meson propagator. 

In general, there is no symmetry between the two arguments of A D in Eq. (6), and 
thus in the rates over the Dalitz plot. A symmetry would be present if, for instance, the 
three-body D decay proceeded only through p-like resonances. We emphasize, however, that 
the product A D (si 2 , s 13 ) A* D (s 13 , s 12 ) in the interference term in Eq. (7) is symmetric under 
the exchange s± 2 <-> Si 3 followed by complex conjugation. This fact is used to simplify the 
analysis. 

The moduli of the D decay amplitude A 12jl3 can be measured from the Dalitz plot of 
the D° — > K s -k~-k + decay. To perform this measurement the flavor of the decaying neutral 
D meson has to be tagged. This can be best achieved by using the charge of the soft pion 
in the decay D* + — > D°ir + . However, the phase #12,13 of the D meson decay amplitude is 
not measurable without further model dependent assumptions. The cosine of the relevant 
phase difference may be measured at a charm factory (see section III). If the three-body 
decay D° — > K s it~it + is assumed to be resonance dominated, the Dalitz plot can be fit 
to a sum of Breit-Wigner functions, determining also the relative phases of the resonant 
amplitudes. This is further discussed in section IV. Here we assume that no charm factory 
data is available and develop the formalism without any model dependent assumptions. 

Using the trigonometric relation cos(a + b) = cos a cos b — sin a sin b, the last term of (7) 
can be written as 

He [A D (s 12 , s 13 ) A* D (si 3 , si 2 ) e^^] = (8) 
^4i2,i3 Mz,i2 [cos(5i2,i3 - £13,12) cos(8 B - l) + sin(5i 2 ,i 3 - #13,12) sin(5 B - 7)] . 

Obviously, to compare with the data, an integration over at least some part of the Dalitz 
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FIG. 1: The partitions of Dalitz plot as discussed in text. The symmetry axis is the dashed line. On the 
axes we have s i2 = m 2 K 7r _ and s i3 = m 2 K 7r+ in GeV 2 . 



plot has to be performed. We therefore partition the Dalitz plot into n bins and define 



J dp A 12 ,iz A 13>12 cos(5i2,i 3 - #13,12), 
/ dp A 12 ,i3 A 13:12 sin(<5i2,i3 - #13,12), 

J i 



Ti — I dp A\ 21Z , 



(9a) 
(9b) 
(9c) 



where the integrals are done over the phase space of the i-th bin. The variables q and Sj 
contain differences of strong phases and are therefore unknowns in the analysis. The variables 
Ti, on the other hand, can be measured from the flavor tagged D decays as discussed above, 
and are assumed to be known inputs into the analysis. 

Due to the symmetry of the interference term, it is convenient to use pairs of bins that 
are placed symmetrically about the 12 <-> 13 line, as shown in Fig. 1. Consider an even, 
n = 2k, number of bins. The k bins lying below the symmetry axis are denoted by index %, 
while the remaining bins are indexed with %. The i-th bin is obtained by mirroring the i-th 
bin over the axis of symmetry. The variables q, Sj of the i-th bin are related to the variables 
of the i-th bin by 

Sj, (10) 

while there is no relation between Tj and Tj. Note that had one used 12 <-> 13 symmetric 
bins centered on the symmetry axis, one would have had Sj = 0. 

Together with the information available from the B + decay, we arrive at a set of Ak 
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equations 



f- = fdt{B- - (K S 7i-7i + ) D K-) = 

Ti + r 2 B T- +2r B [cos(5 B - t)c* + sin (^ - 7)«i], 
ff = [df(B- - (7^7^+)^-) = 

J i 

Ti + r 2 B Ti +2r B [cos(5 B - 7)q - sin(5 B - 7)5;], 
ft = jdt{B+ -> (^ 5 7r-7T+) D K + ) = 

7j + r|Tj +2r B [cos(5 B + 7)^ - sin(5 B + 7)^], 

f-+ EE jdt (B + -> (^ 5 7T-7T + ) D K + ) 



(11a) 



(lib) 



(lie) 



(lid) 



^ + r|T- +2r B [cos(5 B + 7 ) Cl + sin(5 B + 7)^]. 

These equations are related to each other through 12 <-> 13 and/or 7^—7 exchanges. All 
in all, there are 2 A; + 3 unknowns in (11), 

q, 5^, r B , 5b, 7, (12) 

so that the 4k relations (11) are solvable for k > 2. In other words, a partition of the D 
meson Dalitz plot to four or more bins allows for the determination of 7 without hadronic 
uncertainties. This is our main result. 

Alternatively to this binning, one can use a partition of the Dalitz plot into k bins which 
are symmetric under 12 <-> 13. For that case, Sj = and the set of the 4k equations (11) 
reduces to 2k relations (the first two and the last two equations in (11) are the same in this 
case). Then, there are just k + 3 unknowns to be solved for, which is possible for k > 3. 
While such binning may be needed due to low statistics, it has several disadvantages, which 
are further discussed below. 

When d = or s« = for all i, some equations become degenerate and 7 cannot be 
extracted. However, due to resonances, we do not expect this to be the case. Degeneracy 
also occurs if 5 B = 0. In this case, 7 can still be extracted if some of the q and/or Sj are 
independently measured, as discussed in the following sections. 

The optimal partition of the Dalitz plot as well as the number of bins is to be determined 
once the analysis will be done. Some of the considerations that enter this choice are as 
follows. First, one would like to have as many small bins as possible, in order that q and 
Si do not average out to small numbers. Second, the bins have to be large enough that 
there are significantly more events than bins. Otherwise there will be more unknowns than 
observables. There are also experimental considerations, such as optimal parameterization 
of backgrounds and reconstruction efficiency 
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III. IMPROVED MEASUREMENT OF a AND Si 

So far, we have used the B decay sample to obtain all the unknowns, including q and 
Si, which are parameters of the charm system. We now discuss ways to make use of high- 
statistics charm decays to improve the measurement of these parameters, or obtain them 
independently. Doing so will reduce the number of unknowns that need to be determined 
from the relatively low-statistics B sample, thereby reducing the error in the measurement 
of 7. 

The first improvement in the measurement is obtained by making use of the large sample 
of tagged D decays, identified in the decay D* + — > D°tt + , at the B factories. So far we only 
assumed that we use this data to determine Tj. In fact, it can also be used to bound the 
unknowns q and s$ defined in (9): 

N, N < [dp A 12A3 A 13A2 < \fT~Ti. (13) 

J i 

This bound will help decrease the error in the determination of 7, with an especially signif- 
icant effect when, due to low statistics in each bin, q and s$ will be determined with large 
errors. 

Next, we show that the q can be independently measured at a charm factory [14-16]. 
This is done by running the machine at the ^(3770) resonance, which decays into a DD 
pair. If one D meson is detected in a CP eigenstate decay mode, it tags the other D as an 
eigenstate of the opposite CP eigenvalue. The amplitude and partial decay width for this 
state to decay into the final state of interest are 

A(D° ± -> K s ( Pl )ir~(p 2 )ir + (p 3 )) = -±= (A D (s 12 , s 13 ) ± A D (s 13 , s 12 )) , (14) 
dT(D± -> K s ( Pl )n-(p 2 )n + (p 3 )) = ^ (A 2 1213 + A 2 1312 ) ± A 12 , 13 A 13jl2 cos(<5 12 , 13 - 5 13A2 )dp. 
where we defined D± = (D° ± D )/ v / 2. With these relations, one readily obtains 



1 

Ci = - 
2 



JdT{D° + ^ K s (pi)7i-(p 2 )7i + (p 3 )) - JdT{D Q _ ^ K s (pi)7r~(p 2 )n + (p 3 )) 



(15) 



As stated above, obtaining this independent measurements reduces the error in the mea- 
surement of 7 by removing k of the 2k + 3 unknowns. 

We can further improve the measurement if we take each bin i and further divide it into 
rii sub-bins, such that the quantities A± 2t i 3 , cos(5i 2 ,i3 — #13,12), and sin(<5i2,i3 — #13, 12) do not 
change significantly within each sub-bin i'. Naively, this statement appears to introduce 
model dependence. In practice, however, the high statistics in the tagged D sample and the 
charm factory -^(3770) sample allow its verification up to a statistical error, which can be 
measured and propagated to the final measurement of 7. 
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Given this condition, Eq. (9a) may be written as 



Ci = 



(16) 




where the i'-th sub-bin is the 12 <-> 13 mirror image of the i'-th sub-bin, Aj/ and 5j/ are 
the values of A 12)13 and £12,13 on sub-bin i', taken to be constant throughout the sub-bin, 
and Ap^ is the area of sub-bin %' . Analogously to Eq. (9c), we have defined the quantities 
Tj/ = Al 213 Apit, which are measured using the tagged D sample. The Q/'s are assumed 
to be measured at the charm factory, applying (15) to the sub-bin i' . Similarly, Eq. (9b) 
becomes 



Eq. (17) removes the k unknowns Sj, and replaces them with the two-fold ambiguity asso- 
ciated with the sign of the square root. Thus, the best approach is to have the signs of s« 
determined by the fit, while constraining their absolute values to satisfy Eq. (17). Doing so 
will reduce the "strain" on the B decay sample, reducing the error on 7. 

Another option for removing the dependence on Sj is to use bins centered symmetrically 
about the 12 <-> 13 line, making s« vanish, as discussed after Eq. (10). In this case, both the 
number of unknowns and the number of observables (bins) is reduced by k. By contrast, 
using Eq. (17) introduces new information from the independent tagged D sample, and is 
therefore preferred. Doing so also preserves the sin (5b — 7) terms in Eq. (11), which helps 
resolve discrete ambiguities (see [7] and section V). 

IV. ASSUMING BREIT-WIGNER DEPENDENCE 

If the functional dependence of both the moduli and the phases of the D° meson decay 
amplitudes Ad(«i2, S13) were known, then the analysis would be simplified. There would be 
only three variables, tb, 5b, and 7, that need to be fit to the reduced partial decay widths in 
Eq. (7). A plausible assumption about their forms, which is also supported by experimental 
data [17-19], is that a significant part of the three-body D° — > K s it~it + decay proceeds 
via resonances. These include decay transitions of the form D° — > K s p° — > K s it~it + or 
D° — > A'*~(892)7r + — > K s it~it + , as well as decays through higher resonances, e.g., / (980), 
/ 2 (1270), or / (1370), inducing p-like transitions, or K*(U30), which induces a iT*(892)-like 



It is important to stress that these assumptions can be tested. By making use of the 
high statistics tagged D sample, one can test that the assumed shapes of the resonances are 
consistent with the data. While the error introduced by using the Breit-Wigner shapes is 



]T ^W^r sin(5, -5 V ) = J2 ±^T t ,T v -cl 



(17) 



transition. 
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theoretical, it is expected to be much smaller than the statistical error in the measurement 
of 7. It will become a problem only when the B sample is large enough to provide a precision 
measurement of 7. By then the tagged D sample will have increased as well, allowing even 
more precise tests of these assumptions, as well as improving the precision of the methods 
presented in section III. 

The decay amplitude can then be fit to a sum of Breit-Wigner functions and a constant 
term. Following the notations of Ref. [20] we write 

A D (s 12 ,s 13 ) = A(D° - K s ( Pl )n~(p 2 )ir + (p 3 )) = 

= a e i5 ° + J2 a r e tSr A r (s 12 , s 13 ), (18) 

r 

where the first term corresponds to the non-resonant term and the second to the resonant 
contributions. The Breit-Wigner function is defined as 

A r {s 12 ,s 13 ) = J M r x BW r , (19) 

where r represent a specific resonance in either the K s (pi)n~ (p 2 ), K s (pi)n + (p 3 ) or 
it~ {p 2 )i{ + {p 3 ) channel. J M. r is the term which accounts for the angular dependence. It 
depends on the spin J of the resonance. For example, °M r = 1 and 1 M. r = —2k\ ■ k 3 . Here 

— * — * 

ki, k 3 are, respectively, the three momenta of one of the particles originating from the res- 
onance and of the remaining particle, as measured in the rest frame of the two resonating 
particles [20]. BW r corresponds to the relativistic Breit-Wigner function and is given by 

BW " {S) = ,-M? + k.r.W < 20) 

where M r is the mass of the r-th resonance and F r (y/s) denotes the mass-dependent width. 
The argument of BW r is s 12 [s 13 , s 23 ] for a K s (pi)n- (p 2 ) [Ks(pi)n + (p 3 ), n~ {p 2 )n + {p 3 )\ 
resonance. One can find detailed expressions for all the functions mentioned above in Ref. 
[20]. 

One of the strong phases 5i in the ansatz (18) can be put to zero, while others are fit to the 
experimental data together with the amplitudes aj. The best option is to fit the Dalitz plot of 
tagged D decays, as was done a decade ago by the ARGUS and E687 collaborations [17, 18] 
and recently by the CLEO collaboration [19]. The obtained functional form of Ad(si 2 , Si 3 ) 
can then be fed to Eq. (7), which is then fit to the Dalitz plot of the B^ — > (Ks'k~'k + )e>K ± 
decay with tb, 5b and 7 left as free parameters. In appendix B we provide a formula for 
the latter case, where only three resonance are included in the analysis. 
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V. DISCUSSIONS 



The observables tf defined in (11) can be used to experimentally look for direct CP 
violation. Explicitly, 

°cp = r * 7 - = 4r -B sin 7 [ci sin 5 B - s» cos 5 B ] , 

°cp = r i~ - r » + = 4r s sin 7 fa sin 5 B + Si cos 5 B ] • (21) 

It is manifest that finite acp requires non vanishing strong and weak phases. The first terms 
in the parenthesis in (21) depends on sin 5b. This is the same dependence as for a two-body 
D decays into CP eigenstates. In the second terms, which depend on cos 5b, the required 
strong phase arises from the D decay amplitudes. Due to the resonances, we expect this 
strong phase to be large. Therefore, it may be that direct CP violation can be established 
in this mode even before the full analysis to measure 7 is conducted. With more data, 7 
can be extracted assuming Breit-Wigner resonances (cf. section IV). Eventually, a model 
independent extraction of 7 can be done (cf. section II and III). 

The above proposed method for the model independent measurement of 7 involves a 
four- fold ambiguity in the extracted value. The set of equations (11) are invariant under 
each of the two discrete transformations 

Pit = {5b -> S B + 7r,7 -> 7 + tt}, P_ = {5 B -> -S B , 7 -7, s i ~ s i}- (22) 

We note that if all the bins used are symmetric under 12 <-> 13, the absence of the sin(5 B — 7) 
terms in Eq. (11) introduces a new ambiguity transformation, P ex = 7 — > 5 B ,5 B — * 7. The 
discrete transformation P n is a symmetry of the amplitude (6) and is thus an irreducible 
uncertainty of the method. It can be lifted if the sign of either cos^b or sin5 B is known. 
The ambiguity due to P_ can be resolved if the sign of sin5 B is known or if the sign of Sj 
can be determined in at least some part of the Dalitz plot. The latter can be done by fitting 
a part of the Dalitz plot to Breit-Wigner functions. We emphasize that only the sign of the 
phase of the resonance amplitude is required, and thus we can safely use a Breit-Wigner 
form for this purpose. 

The r B suppression present in the scheme outlined above can be somewhat lifted if the 
cascade decay B~ — > DX~ — > (K S Tr~n + ) D X~ is used [6, 11]. Here X~ is a multibody 
hadronic state with an odd number of kaons (examples of such modes are K~it~it + , K~n 
and i^57r _ 7T°). Unlike the B~ — > D°K~ decay, these modes have color-allowed contributions. 
This lifts the color suppression in r B , while the mild suppression due to the CKM matrix 
elements remains. The major difference compared with the case of the two-body B~ decay 
is that now r B and 5 B are functions of the B — > DX~ decay phase space. Therefore, the 
experimental analysis has to deal with two Dalitz plots, one describing B — > DX~ and the 
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other describing the D — > Ksn~n + decay. In appendix C the necessary formalism that 
applies to this case is outlined. Note that the above mentioned treatment for multibody B 
decays also applies to quasi two-body B decays involving a resonance, such as B — > DK* . 

In addition to using different B modes, statistics may be increased by employing various 
D decay modes as well. An interesting possibility is the Cabibbo allowed D — > K s tt~it + it 
decay. It comes with an even larger branching ratio than the D — > Ksii~i{ + decay. In 
addition, it has many intermediate resonances contributing to the greatly varying decay 
amplitude, which is what is needed for the extraction of 7. The disadvantages of this mode 
are the low reconstruction efficiency of the 7r°, as well as the binning difficulties introduced 
by the higher dimensionality of the four-body phase space. The formalism of section II 
applies to this mode as well, but now the partition of the four-body phase space is meant 
in Eq. (11). In the equivalent of (5), this mode has an extra minus sign, since we have 
introduced a new CP-odd state, the 7r°. The final set of equations is then obtained from 
(11) by replacing r B — > — r B . The Cabibbo allowed mode D — > K~K + K S may also be used 
for the extraction of 7, as can the Cabibbo suppressed decays to K~K + ir°, 7T~7r + 7r°, and 
K s K + ir~. One can also use (almost) flavor eigenstate decay modes, such as D — > K~it + it 
and D — * K~it + it~it + [5]. Here, the important interference is between the Cabibbo allowed 
D decay and the doubly Cabibbo suppressed D decay. 

While we concentrated on charged B decays, the Dalitz plot analysis presented here can 
also be applied to self-tagging decays of neutral B mesons [8]. It is also straightforward to 
apply it to cases where time dependent CP asymmetries are measured [2]. 

The sensitivity to 7 is roughly proportional to the smaller of the two interfering ampli- 
tudes. Assuming that the only two small parameters are r B and A, our method is sensitive 
to 7 at 0(r B ). However, the method is sensitive to 7 only in parts of the Dalitz plot. The 
highest sensitivity is in regions with two or more overlapping resonances. The sensitivity of 
the proposed method is therefore of order 0(r B £), where £ 2 is the fraction of events which 
are in the interesting region of the Dalitz plot. 

A crucial point of our method is that it uses interference between two Cabibbo allowed D 
decay amplitudes. This is against the common intuition, which suggests that we must have a 
A 2 suppression for such interference to take place, as we need a final state that is common to 
both D and D. Specifically, one typically requires one Cabibbo allowed decay and another 
that is doubly Cabibbo suppressed, or two decays that are singly Cabibbo suppressed. To 
overcome this preconception, our method makes use of K° — ~K° mixing (which is also the 
case for the two-body D — > K s ir° decay), plus the existence of overlapping resonances, which 
are obtained by Cabibbo-allowed D° and D° decays. In addition, it is important that the 
hadronic three-body D meson decays have a widely changing amplitude over the Dalitz plot, 
which is ensured by the presence of resonances in this energy region. If the strong phases 
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#12,13 and the moduli A 12> i3 in Eq. (9) were (almost) constant across the available phase 
space, the extraction of 7 from Eqs. (11) would not be possible. 

Before concluding, we mention that quasi two-body D decays where one of the particles 
is a resonance, such as D — > K* + ir~ and D — > K + p~ [4], were proposed for use in measuring 
7. But in fact, using such decays requires a Dalitz plot analysis (see e.g. [10, 12]). What 
we showed here is that one can actually use the whole Dalitz plot to carry out the analysis 
and does not need to single out contributions of one particular resonance. Moreover, we 
showed that the assumption about the shapes of the resonances can be avoided, essentially 
with currently available data-sets. 

In conclusion, we have shown that the angle 7 can be determined from the cascade decays 
— > K ± (K s tt~tt + ) d . The reason for the applicability of the proposed method lies in the 
presence of resonances in the three-body D meson decays that provide a necessary variation 
of both the phase and the magnitude of the decay amplitude across the phase space. The 
fact that no Cabibbo suppressed D decay amplitudes are used in the analysis is another 
advantage of the method. However, it does involve a Dalitz plot analysis with possibly 
only parts of the Dalitz plot being practically useful for the extraction of 7. In reality, 
many methods have to be combined in order to achieve the required statistics for a precise 
determination of 7 [7]. 
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APPENDIX A: THE EFFECT OF D-D MIXING 



In this section we focus on the contributions introduced by the fact that the flavor states 
|-D°), |Z?°) and the mass eigenstates \D h ,l) — Pd\D°) ±q D \D°) do not coincide. This effect 
was studied in the general case in Ref. [21]. Here we apply their formalism to our case. 

Following Ref. [21] we introduce the rephasing-invariant parameter xi 

1 + AD^fi^B-^D 
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where 

Ad-»/ = = ^ = rBe-*^^^, (A2) 

and we use the definitions of Eqs. (3) and (4) and allow for new physics effects in qn/PD — 
e l20D . (In the phase convention where the D decay amplitudes are real, the phase 6d is 
negligible in the Standard Model). In our case, the final state / equals K S 7i~7r + , which 
leads to 

*d-,k s { Pi )-k-(p 2 )k+(p z ) - e — t- — — - n D {si2,, si3je '. (Adj 

Ad(S12, Si 3 ) 

Once D — D mixing is taken into account in the analysis, the expression for the partial decay 
width (7) is multiplied by the correction term [21] 

l-ne( X i)yD+Im(xi)x D , (A4) 

where we have expanded the correction term to first order in the small parameters 

Am AT . . 

xd = — , y D = ^r, (A5) 

where Am and Ar are the mass and decay width differences in the D — D system, and T is 
the D° decay width. The values of xd and yo are constrained by present measurements to 
be in the percent range, yo = (1-0 ± 0.7)% [22] and \x\ < 2.8% [23] (assuming small strong 
phases) . 

The ratio of magnitudes, Rd(su,, S13), depends on the position in the Dalitz plot, and 
can vary widely. Our method is useful for the model independent extraction of 7 only in 
the region where Rd is of order one. We therefore distinguish three limiting cases 

• Rd 3> 1 r B , for which TZe(xi),^ m {Xi) ~ 0(l/r B ) and therefore the corrections in 
(A4) can be of order 10%. However, this is the region of Dalitz plot where our method 
is mostly not sensitive to 7 and therefore the induced corrections due to D — D mixing 
do not translate into an error on the extracted 7. 

• Rd ~ 1 > r B , for which 1Ze(xi) ,1m(xi) ~ 0(1) and therefore the corrections in (A4) 
are at the percent level. This is the value of Rd for which our method is most sensitive 
to 7. 

• 1 ^> r B ~ Rd, for which TZe(xi) ^ m (Xi) ~ 0(r B ,R B ) and therefore the corrections 
in (A4) are very small. 

In conclusion, we expect errors of at most a few percent due to neglecting D — D mixing in 
our method. In principle, even these errors can be taken into account [16, 21, 24]. 
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APPENDIX B: A FIT TO BREIT-WIGNER FUNCTIONS: AN ILLUSTRATION 
FOR THREE RESONANCES 

In this appendix we provide the formulae for the fit of D meson decay amplitude to a 
sum of three Breit-Wigner functions describing i^* ± (892) and p° resonances. We write Eq. 
(18) explicitly as 

A D (s 12 ,s 13 ) = A(D° - K s ( Pl )7r~(p 2 )ir + (p 3 )) = 

(Bl) 

= a p A p o(s 23 ) + a K *e lSF A K *(s 12 ) + a K *r D e lSD A K *(s 13 ), 

where 5 F (5 D ) is the strong phase of the Cabibbo favored (doubly Cabibbo suppressed) 
D° -> K*-n+ (D° -> K*+n-) decay with respect to the decay D° -> K s p°. We further 
introduced 

a p oc A(D° - p°K s ) = A(D° - p°K s ), 
a K *e i5F oc A(D° -> IT"^*) = A(D° -> TT+tt-), 
a^r D e lfe oc A(D° -> TT+tt-) = A(D° -> iT-7r + ). (B2) 

The Breit-Wigner functions „4 r are defined in (19), where we write in (Bl) only the s a b 
dependence of the BW r part, given in (20). The first index of s a b is understood to denote 
also the particle appearing in the expression for 1 Ai r (19). Exchanging a <-> b corresponds 
to 1 A4 r <-> — \M r , in particular *4. p o(s 2 3) = _ A p o(s 32 ). In the above we assumed that there 
is no CP violation in the D decays amplitudes. Note that there are two small parameters 

r B ~0.1-0.2, r D ~A 2 ~0.05. (B3) 

We then obtain (cf. (6)) 

A{B~ - (K s (pi)n~(p 2 )n + (p 3 )) D K-) = (B4) 
A B V D x { (a p A p o(s 23 ) +a K * [e iSp A K ,{s 12 ) + r D e l&D A K *{s lz )]) + 

r B e t(&B -^ {a p A p o(s 32 ) + a K * [e^A K '{s 13 ) + r D e idD A K *(s 12 )] ) }. 

The corresponding expressions for B + decays are obtained by changing 7 — > —7 and 
Tr-(p 2 )-K + (p 3 ) -> 7r + (j9 2 )vr"(p 3 ). 
We further define 

5- = aig[A K *(s 12 )], 5 + = Bxg[A K *(si 3 )}, S = &rg[A p o(s 23 )]. (B5) 
where the dependence of 5±^ on the position in the Dalitz plot is implicitly assumed. The 
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reduced differential decay rate is then 

dT(B- -> (K s n'7i + ) D K-) oc 

a 2 p \A p o(s 23 )\ 2 (1 - 2r B cos(5 B - 7) + r|) + 

(s 12 ) | 2 (1 + 2r B r D cos(5f D - 7) + (r B r D ) 2 ) + 
a^, l^x-. (si 3 ) | 2 (rf> + 2r B r D cos(5 BF - 7) + r|) + 

2a p a^»|^o(s23)^4x*(si3)| x 

|r D cos + - r| cos 5q + - r B r D cos(5f h - 7) + r B cos(5<f F+ + 7) | + 
2a p a K *\A p o(s 2Z )A K *{s 12 )\ x 

I cos 8q~ - r B cos(5^ Q - 7) + r B r D cos(5 B d_ + 7) - r|r D cos j + 

2a 2 ^|„4x*(si2)-4x*(si3)| x 

|r D cos 5£+ + r B cos(5 B+ + 7) + r B r 2 D cos(5+_ - 7) + r\r D cos 5^+ } , (B6) 

where the notation of the strong phases is such that the lower (upper) indices indicate phases 
appearing with a plus (minus) sign. For example, 

Sot = 5 D + 5--5 F - S+. (B7) 

a p , clk* and t\d are assumed to be known and thus there are five unknowns to fit, namely 

r B , S D , S F , S B , 7- (B8) 

Using both B~ and B + decays, there is enough information to determine them all. This is 
true even if one neglects terms that scale as r B and even if r B = 0. This indicates that the 
method does not rely on doubly Cabibbo suppressed decays of the D, and that it is sensitive 
to 7 in terms of order r B , rather than r B (See discussion in [10]). Moreover, even if some 
or all of the strong phases that arise from two-body decays, namely, S B , $d, and Sf, vanish, 
there is still enough information to determine 7. 

APPENDIX C: MULTIBODY B DECAY 

We consider the cascade decay B~ — > DX~ — > (Ks7t~7t + )dX~ . Let us assume that the 
phase space of the first decay, B~ — > DX~, is partitioned into m bins that we label by the 
index j, and the phase space of the D meson decay is partitioned into n = 2k bins labeled 
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by % and % as in section II. Instead of Eqs. (11) we now have the set of Ak x m equations 
f". = / dT(B- -> (K s n-n+) D X;) = 

Ji.i 



(Cla) 



Tt + RfTj + cos-fiacf + Si sf) + sin 7 (Qsf - s iC f), 



rr. 



(cib) 

T T + RfTi + cos 7 (qcf - Si sf) + sin 7 (qsf + Si cf ), 



(Clc) 

T T + ^f^ + cos 7 (cicf - s iS f) - sin 7 ( Cl sf + Si cf ), 



f- + = / dT{B + - (7^^+)^+) = 

J Aj (cid) 

Tj + RfTj + cos 7 (qcf + Sjsf ) - sin 7 (Qsf - s^cf ), 

where the integration is over the phase space of the j-th bin in the B decay and the phase 
space of the i-th bin in the D decay. The j-th bin of the B + decay phase space is obtained 
from the j-th bin of the B~ decay by CP conjugation. We also used 



sf = / 2r B sm5 B , 



3 



cf = I 2r B cos5 B , 



-I 



"3 



■I 



Rf = Jfi, (C2) 

where r B and 5 B are functions of the position in the B decay phase space. From the set 
of Ak x m equations (CI), one has to determine 2k + 3m + 1 unknowns q, Sj, cf , sf , Rj , 
and 7 . With a partition of the D decay phase space into 2k > A bins and with a partition 
of the B decay phase space into m > 1 bins, one has enough relations to determine all the 
unknowns, including the angle 7 . This is true even for constant S B and r B , in which case the 
above equations fall into Ak sets of m equivalent relations, i.e. the set of Ak x m equations 
is reduced to the set of Ak independent relations (11). 

Finally, we note that the above equations can be used to determine 7 also for two-body 
D decays [6]. 
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